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Abstract. In this paper we report a new method for determining the thickness of the non-edge-on disk galaxies. 

Our method is based on the observation comparing with the theoretical researches of the distribution of the 
vertical velocity dispersion, which is obtained from the solutions of three dimensional Poisson’s equation and the 
galactic dynamical equation. This method also allows us to investigate the mass-to-light ratio of the disk. As 
examples, the thickness and mass-to-light ratio of two disk galaxies, NGC 1566 and NGC 5247, which have been 
extensively studied in spectroscopy, have been calculated and the results are presented. 
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1. Introduction 

The thickness of disk galaxies is a very important pa¬ 
rameter for understanding these objects, although all of 
which are very thin comparing with the length of the 
disk. Sancisi and Allen (1979) estimated the thickness of 
the edge-on Sb galaxy NGC 891, on the basis of observa¬ 
tion of neutral hydrogen. Van der Kruit & Searle (1981a) 
proposed a model for light distribution in the disks of 
edge-on spiral galaxies, assuming that a galaxy has a lo¬ 
cally isothermal, self-gravitating, truncated and exponen¬ 
tial disk. The model has the feature of being isothermal 
in z-direction at all radii with a scale parameter of z o and 
has an exponential dependence of surface brightness on r 
with scale length of ra- The space luminosity of this model 
can be described by 

L(r , z) = Lqe~ r ^ rd sec h 2 (z/zff) (1) 

With this model, van der Kruit & Searle (1981a, 1981b, 
1982a, 1982b, named KS hereafter) determined and zq 
for several edge-on galaxies without an appreciable bulge. 
Unfortunately, this method may not be suitable for gen¬ 
eral disk galaxies, e.g. face-on galaxies. 

Peng (1988) put forward a method to measure the 
thickness of face-on galaxies on the basis of asymptotic 
formula of the disturbed gravitational potential. And a 
revised method based on the exact integral expression for 
disturbed gravitational potential has been presented by 
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Zhao et al. (2004) to estimate the thickness of face-on 
disk galaxies. But the results obtained by these two cor¬ 
relative methods should be checked with the results gotten 
by other independent ways. 

In this paper, we present the solution of the Jeans 
equation along ^-direction first and then describe how 
to use this solution to determine the thickness of face- 
on disk galaxies. In Section 3, we show the applications 
of our method to two near face-on galaxies, NGC 1566 
and NGC 5247, which have been extensively studied in 
spectroscopy. The results are also presented. And a brief 
discussion is given in the last section of this paper. 


2. Method 


2.1. The distribution of the vertical velocity dispersion 


For a stable, axisymmetric galaxy, the Jeans equation 
along the ^-direction is 


dz 


(p(v?)) 



( 2 ) 


where ( V?) is square of the ^-direction velocity dispersion. 

As the first of the fundamental assumptions, we accept 
Parenago’s density distribution law along the ^-direction 
for a finite thickness galaxy, 


p(r^)=p(r,0)e-^l = |a(r)e-“N, 


( 3 ) 


where a is a parameter the reciprocal of which is half of 
the effective thickness, H, of a galaxy, a may be taken as 
a constant, for the scale height is basically independent of 
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the radius at least for the thick disk of our Galaxy (e.g. 
see Freeman, 1987). cr(r) is its surface density, de Grijs & 
van der Kruit (1996) have also confirmed that the simple 
exponential fit of the vertical stellar light distribution is a 
good approximation. Hence, the potential caused by this 
density can be found by solving the Poisson’s equation of 
the potential 

V 2 </?(r, z) = 27rGaa(r)e~ a ^ z i (4) 


The rigorous solution of the potential has been given 
by Peng et al. (1978) 


<p(r, z) = -7rGa 

'/o°° - ae-0 |z| )J°(/3r)s(/?)d/? 


where 
s(/3) = 



rJo((3r)cr(r)dr. 


( 6 ) 


And Jo((3r) is the well known Bessel function of order 0. 

Then using equation (5), we integrate equation (2) 
from 2 to oo (z > 0) at both sides of the equal sign re¬ 
spectively, it follows (Huang, Huang & Peng 1979) 


(p( V z))z> o = -27rGa 2 p(r, 0) 


f°° H i 

J0 f3 2 -a 2 a+/3 c 


(ct+0)z 


1 p -2 az 
2a c 


J 0 (f3r)s((3)dp 


(7) 


Hence, on the galactic plane (z = 0 ), {V 2 ) could be 
written as 


<K 2 >( 2= 0 ) = nGa / 0 °° T -l w J 0 (f3r)s(P)d/3 
where 

P(r,a) = ^ 7 y/ 0 °° (3((a, f3)Jo(f3r)s((3)d(3, 
C(a,/3) = 1 - (1+/3 1 / a )a • 


( 8 ) 


(9) 


2.2. The thickness parameter a 

By assuming that a galaxy has an infinitesimally thin disk, 
Freeman (1970) studied 36 spiral and SO galaxies with sur¬ 
face photometry and showed that the radial light distribu¬ 
tion in the disks of spiral galaxies can be described by an 
exponentially decreasing surface-brightness with increas¬ 
ing galactocentric radius. Along the vertical direction, de 
Grijs & van der Kruit (1996) have shown that a simple 
exponential fits turned out to be good approximations of 
the stellar light distribution. Thus, the model to account 
for the space-luminosity can be described by 

L(r,z) = L 0 e- r/rd e~ alzl , (10) 

where r is the position along the major axis, Lq the central 
space luminosity in the plane of the galaxy, r<i the scale 
length of the disk. From this we calculate the face-on dis¬ 
tribution of surface-brightness 

I(r) = I 0 e- r / rd , (11) 


where Iq is the central surface brightness of the disk. 

Making use of the observed constant color index (Van 
der Kruit & Searle 1981b, 1982a), Van der Kruit & 
Freeman (1986) pointed out that the mass-to-light ratio of 
the old disk , Y* = ( M/L) 0 id disk, is approximately con¬ 
stant with the radius. This means we’ll take the projected 
surface density as 

a(r) = T*/(r) = TJ 0 e~ r/rd = a 0 e~ r ^, (12) 


where do is the surface density at r = 0. 

Therefore, equation (6) could be reduced by substitut¬ 
ing a(r) with equation (12) to 


r°° r 2 

m-l r Jo {(dr) a (r) dr = T.J 0 ^ + r ^ 2)3/2 • ( 13 ) 

And the expression of (V 2 ) may be rewritten as 
<K 2 > = ^T*/(r)[l-P(r)] 

(14) 

= 2fTJ 0 [e- r / rd -P 1 {r)] 

where 


Pi(r) = P(r)e~ r / rd 

= Jo (T+ 7 p 2 j 372 ^[l - (i+p/a) 2 \J°(P r )dP 


(15) 


Hence, if we know the /o, Td from photometric study 
and (V 2 ) from spectroscopic observation, we can calcu¬ 
late the thickness parameter a through equation (14) and 
therefore obtain the thickness of galaxies. 


3. Application 

We found that there have existed extensively spectroscopic 
study on two nearly face-on galaxies NGC 1556 and NGC 
5247. Therefore, we choose these two galaxies as examples 
to illustrate how our method works. The general informa¬ 
tion of these two galaxies are summarized in Table 1. 

3.1. Photometric decomposition 

Although there also existed photometric results for NGC 
1566 (de Vaucouleurs 1973; Bottema 1992) and 5247 (van 
der Kruit & Freeman 1986), we can not use these data 
directly since the photometric decomposition is very im¬ 
portant for our purpose. The photometric decomposition 
can disentangle the contributions to the total luminosity 
of the bulge and disk, and this will allow the identification 
of the disk region. Moreover, because the infrared images 
would suffer much less extinctions than optical ones do, 
we use the 2MASS H-band images to do photometric mea¬ 
surements. 

To obtain photometric parameters, and /jq, we 
use the two dimensional image-fitting algorithm GALFIT 
(Peng et al. 2002) designed to extract structural param¬ 
eters directly from the galaxy images. GALFIT assumes 
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Table 1 . Parameters of NGC 1566 and NGC 5247. 


NGC 

Inclination 

Classification 

Distance 0 (Mpc) 

M h (mag) 

ref 

1566 

28° ± 5° 

Sc(s)I, Seyfert 1 

17.4 

-24.0 

(1), (2), (4) 

5247 

20° ± 3° 

Sc(s)I-II 

16.0 

-23.2 

(1), (3), (4) 


a Obtained based on Ho=75 km s 1 Mpc 1 . 

References (1) Sandage & Tanmiann 1981; (2) Bottema 1992; (3) van der Kruit & Freeman 1986; (4) Jarrett et al. 2003. 



Fig. 1. For both galaxies we show the original images (left), the smooth and symmetric model images fitted by 
GALFIT ( middle ) and the residual images (right). For NGC 1566 (top), the size of each field is 2'.4 x 2'.9; for NGC 
5247 (bottom), the size is 6'.7 x 6'.7 


a two-dimensional model profile for the galaxy. The func¬ 
tional forms of the models we choose to fit include combi¬ 
nations of a Gaussian, a Sersic r 1 '" law and an exponen¬ 
tial disk profile. For the case of a Seyfert 1 galaxy, NGC 
1566, we use all three functions to fit, while we use the 
last two functions for NGC 5247. We fit the following: the 
(x, y) position of the center, M to t (the total magnitude 
of the component), r e (the effective radius), n (the Sersic 


index), q (the axis ratio), r<i (the scale length of the ex¬ 
ponential disk), the major position angle (PA), and c (the 
diskiness/boxiness index). The detailed parameters of all 
fitted components are listed in Table 2. The original im¬ 
ages, the smooth and symmetric model images fitted by 
GALFIT and the residual images are shown in the left, 
middle and right panels in Figure 1, respectively. 
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Table 2. Bulge-disk decompositions for NGC 1566 and NGC 5247. The fitting functions are a Gaussian (for the 
nucleus), Sersic r x ! n law (for the bulge) and an exponent (for the disk). 


Component 

Mtot 

(mag) 

FWHM/r e /r/ 

(kpc) 

Me/hO 6 

(mag arcsec -2 ) 

PA 

(°) 

n 

q 

C 

xl 

NGC 1566 

Nucleus 

ll.70i0.03 

0.32A0.01 


88.02i2.36 


0.89i0.01 

0.26i0.11 


Bulge 

9.25i0.02 

0.92i0.03 

17.05±0.06 

12.21i0.48 

1.65i0.03 

0.79i0.00 

-0.28i0.01 


Disk 

7.65A0.00 

3.29A0.04 

17.20±0.04 

-15.62i0.32 


0.73i0.00 

-0.29i0.02 

0.310 

NGC 5247 

Bulge 

10.69A0.01 

0.67A0.02 

18.06±0.04 

42.20il.44 

1.24i0.01 

0.87i0.00 

0.01i0.04 


Disk 

7.84A0.00 

4.78i0.03 

18.60±0.03 

-18.72i0.70 


0.80i0.00 

-0.08i0.02 

0.130 


a For the nucleus it is FWHM, for the bulge it is r e , and for the disk it is 
b For the bulge it is fj, e , and for the disk it is po- 
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Fig. 2. 2MASS H-band radial surface brightness profiles. Symbols show the radial profiles obtained from ELLIPSE, 
the dotted line shows the Seyfert core component fitted to the Gaussian function, the dashed lines show the bulge 
components fitted to the Sersic r 1 /™ law, and the dashed-dotted lines show the disk components fitted to the exponential 
law. The thick solid lines are the superposition of all the components. The lower portion of each panel shows the 
difference between the data (ELLIPSE) and the fits (GALFIT). 


In order to check whether or not the fitted re¬ 
sults are reliable, we examine the residual images ob¬ 
tained by subtracting the fitted models. In both cases 
we found no obvious evidence (except for the expectant 
spiral arms) for components beyond the fitted models. 
Furthermore, we also measure surface brightness profiles 


using the I RAF 1 / ELLIPSE algorithm. In Figure 2 we 
show the radial profiles obtained from the isophote fit- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which is operated by Associated of Universities 
for Research in Astronomy, Inc., under cooperative agreement 
with the National Science Foundation. 
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ting by ELLIPSE and the profiles of the chosen model 
for each galaxy. The differences between the observa¬ 
tional data ( ELLIPSE ) and the models (GALFIT) range 
mostly over ±0.2 mag arcsec -2 , except for the regions dis¬ 
torted by the spiral arms at large radii. So we conclude 
that our global fits are reliable. 

3.2. Fitting the observed dispersion data 

The vertical velocity dispersions of NGC 1566 were got¬ 
ten from Table 4, column 4 in Bottema (1992); NGC 
5247 from Table 1, column 4 in van der Kruit & Freeman 
(1986). As these two galaxies are not exactly face-on (see 
Table 1), this will affect the observational data at two as¬ 
pects: (1) the observed velocity dispersion is not exactly 
equal to the vertical velocity dispersion. The observed dis¬ 
persion is approximately 4% higher than the real vertical 
dispersion for NGC 1566 (Bottema 1992), ~5% larger for 
NGC 5247 (van der Kruit & Freeman 1986). (2) the ob¬ 
served surface brightness only amounts to about cosi times 
the face-on brightness. Before using the photometric and 
spectroscopic data for fitting, we have corrected the affec¬ 
tions of the inclination. We’ll present a further discussion 
about how the inclination affects the fitted results in the 
following. 

To obtain the thickness parameter a, we use a least- 
squares method to fit the observational data shown in 
Figure 3, and 4, for NGC 1566 and NGC 5247, respec¬ 
tively, and only the points with radius larger than ~ 
2 — 3r e , where the bulge’s contribution to the velocity 
dispersion may be neglect according to Terzic & Graham 
(2005) (see Figure 11 and 12 in their paper). In the pro¬ 
cess of fitting, we first adopt Y*^ (the mass-to-light ra¬ 
tio for H-band) as an empirical value, following Devereux 
et al. (1987), Oliva et al. (1995) and Giovanardi & Hunt 
(1996), and then use equation (14) to derive the thickness 
parameter a and y 2 . Whereafter, we adopt this fitted a 
to derive the mass-to-light ratio Y*^. Such an iteration 
are repeated until Ay 2 /y 2 (Ay 2 is the difference between 
the xt given by two successive fitting) is less than several 
times 10~ 4 . The parameters needed in our fitting process 
and the fitted results are listed in Table 3. The best fitted 
results are also shown in Figure 3 and 4 with dotted lines. 

The values of y 2 given by our least-squares fitting for 
NGC 1566 and NGC 5247 are 0.48 and 0.22, respectively 
(see Table 3). These values of y 2 are somewhat smaller 
than statistically acceptable, presumably because of the 
large measurement errors of the observed data and the 
small samples used for the fittings. However, the best- 
fitted results with our model are very exciting, as seen 
from Fig 3 and 4. They show that our method may be 
realistic and could be used universally to measure the 
thicknesses of non-edge-on disk galaxies. What’s more, 
our method can derive the mass-to-light ratio simultane¬ 
ously. But the value gotten by this way are overestimated 
since the mass contribution of the dark matter halo is 
neglected. However, this is a reasonable assumption. The 


kinematic data cover a radial range where luminous mat¬ 
ter can be assumed to dominate the galaxy dynamics. As 
shown in Table 3, the derived mass-to-light ratio for these 
two galaxies are well consistent with the empirical values 
(Devereux et al. 1987; Oliva et al. 1995; Giovanardi & 
Hunt 1996). This also suggests that our method is reason¬ 
able. 

When this method is applied to other measurements, 
we maybe adopt an easier way that we only need to know 
the vertical velocity dispersion at some exact radius (must 
be larger than ~ 2r e ), (V^) 1 / 2 , and the photometric 
results, e.g. the scale length of the disk rd and the central 
extrapolated surface brightness Iq. Then, equation (14) 
can be rewritten as 

(V* ri )a + nGrj o P 1 (r 1 ,a)-nGrj o e- ri / rd =0 (16) 

The thickness parameter, a, is the only unknown param¬ 
eter in this equation if we have assumed the mass-to- 
light ratio. Alternatively, we can use the surface density- 
weighted average velocity dispersion, (V?), to estimate the 
thickness parameter a through the following equation, 

r* 2^ _ p( r) w r 

(V?) = J ° a Jr w • (17) 

Jo v{r)rdr 

When equation (17) is used for measurements, equation 
(12) would be substituted for er(r). 

4. Discussion 

So far, we have proposed an method to determine the 
thickness of the non-edge on disk galaxies H , based on the 
photometric result, scale length of disk (r^) and extrap¬ 
olated central surface brightness (Jo); and spectroscopic 
result, the vertical stellar velocity dispersion (( V) 2 )) ■ The 
fundamental assumption required to derive the conclu¬ 
sions presented here is that the disk mass profile is expo¬ 
nential, both along the vertical direction (i.e., ^-direction, 
see equation (3)), which is different from the isothermal 
model used by KS, and along the radial direction (i.e., re¬ 
direction, see equation (1)). It has been shown by de Grijs 
& van der Kruit (1996) that the best fitting vertical model 
is more peaked than expected for an isothermal sheet dis¬ 
tribution and the simple exponential fits turned out to 
be good approximations of the vertical surface brightness 
profile. Peletier & de Grijs (1997) also find that the K- 
band light (often assumed to trace mass) indeed has a very 
peaky, almost exponential, vertical distribution in edge-on 
spirals. By contrast, Gould, Bahcall & Flynn (1996) find 
that the vertical distribution of Galactic M dwarfs (which 
probably trace the stellar disk) is somewhat less peaky. 

Because our method is intended to work for a pure 
disk galaxy, the photometric decomposition is another 
important ingredient affecting the fitted result. As de¬ 
scribed in Section 3.1, photometric decomposition allows 
us to identify the real disk region and the range that 
the bulge/nucleus (if the galaxy has) affects. We also use 
the data just outside the effective radius r e to fit a for 
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Table 3. Photometric parameters and least-square fitted results for NGC 1566 and NGC 5247. 


NGC 

Io 

T» 


a 

H a 

Band 


( 1/0 pc -2 ) 

(Mq/Lq) 

xl 

o 

i 

(kpc) 


1566 

1193.6±29.7 

0.49±0.03 

0.48 

1.75A0.56 

1.14±0.36 

H 

5247 

328.7±4.5 

0.91A0.03 

0.22 

1.32±0.52 

1.52±0.60 

H 


a Calculated according to H = 2/a. 
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Fig. 3. Observed (filled triangles) and fitted (dotted line) distribution of vertical stellar velocity dispersion of the disk 
in NGC 1566. The observed data are from Bottema (1992), but the data which belong to the Seyfert core and the 
bulge are not included in the figure, and not taken into account for the determination of the thickness parameter a. 
The line is our least-squares fitted results with equation (14) in the text. 


NGC 1566, and we obtain a much bigger y 2 of 1.8 and 
a smaller a of 0.96 kpc -1 . For NGC 5247, we use the 
data out to the center since there are no data just near 
the effective radius. This also gives a bigger y 2 of 0.35 
and a smaller a of 1.0 kpc -1 . This suggests that we need 
to select the radial range according to the total mass of 
the bulge/nucleus. The inclination is another ingredient 
affecting the fitting results, particularly true for the mass- 
to-light ratio. Applied corrections for different inclinations 
will result in different mass-to-light ratio, but almost the 
same thickness parameter. For the same velocity disper¬ 
sion data, a smaller inclination gives a larger mass-to-light 
ratio. 

The parameter H should be regarded as about twice 
the scale height of the best fit of an exponential to the 
actual distribution. We’ve compared our some relevant re¬ 


sults, e.g. rd/H (the scale length to the effective height ra¬ 
tio), H (the effective height) of these three galaxies with 
those of eight edge-on spirals of KS, which are directly 
measured by fitting the photometric data. As shown in 
Table 4, our results are very consistent with KS. This 
make us believing that our new method for determining 
the thickness of non-edge-on disk galaxies is reasonable 
and feasible. In Figure 5 we plot the thickness versus the 
Hubble type. The continuous line is the fitted result for 
the crosses, which represent the median value of the thick¬ 
ness for each type. The trend that the thickness of disks 
decreases along the Hubble sequence is also found by Ma 
(2002) and Zhao, Peng & Wang (2004). However, it needs 
further study on this correlation because the sample used 
here is rather small. 





Zhao et al.: New method for determining the thickness of disk galaxies 


7 





A 

Cd 

> 

V 


80 

70 

60 

50 

40 

30 

20 

10 



Radius (kpc) 


Fig. 4. Observed (filled triangles) distribution and our fitted results (dotted line) with H-band photometric data for 
the disk of NGC 5247. The observations are obtained by van der Kruit & Freeman (1986). But the data belonging to 
the bulge are excluded from this figure and not taken into account for our fit. 

Table 4. Parameters of the near face-on galaxies from our study compare with those of edge-on galaxies from the 
study by van der Kruit & Searle. 


NGC 

RC3 Type 

M H a 

r d b (kpc) 

H c (kpc) 

r d /H 

ref 

1566 

Sbc 

-24.0 

3.3 

1.1 

3.0 

our study 

5247 

Sbc 

-23.2 

4.8 

1.5 

3.2 

our study 

891 

Sb 

-23.3 

4.9 

0.99 

5.0 

(1) 

4013 

Sb 

-22.4 

3.4 

1.1 

3.1 

(2) 

4217 

Sb 

-22.9 

3.5 

1.7 

2.1 

(2) 

4244 

Scd 

-20.7 

2.6 

0.58 

4.5 

(3) 

4565 

Sb 

-23.7 

5.5 

0.79 

7.0 

(3) 

5023 

Scd 

-19.5 

2.0 

0.46 

4.3 

(2) 

5907 

Sc 

-23.1 

5.7 

0.83 

6.9 

(2) 

7814 

Sab 

-23.5 

8.4 

2.0 

4.2 

(4) 


a The magnitudes are all gotten from Jarrett et al. 2003; 
b r d in I<S is represented by h; 

c H is represented by zo, which is obtained by fitting the light distribution with z) = /r(0, 0)(R/h)Ki(R/h.)sech 2 (z/zo), 
K\ is the modified Bessel function, in KS. 

References (1) van der Kruit & Searle 1981b; (2) van der Kruit & Searle 1982a; (3) van der Kruit & Searle 1981a; (4)van der 

Kruit & Searle 1982b. 


The method presented in this paper is rather difficult 
to be used for measurements of large galactic samples 
because the distributions of the vertical stellar velocity 
dispersions are inconvenient to be obtained in practice. 
Fortunately, the most essential purpose of the present pa¬ 
per is to test the availability of the method proposed by 


Zhao et al. (2004), which is on the basis of the density wave 
theory. Hence, we measured the thickness of the disks of 
NGC 1566 and NGC 5247 used the method put forward by 
Zhao et al (2004) and obtained 0.6 ± 0.4 kpc and 0.7 ± 0.5 
kpc, for NGC 1566 and NGC 5247, respectively. These 
results are somewhat smaller than the results obtained 
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Fig. 5. Thickness plotted versus Hubble type(l, 2, 3, 4, 5 and 6 for Sa, Sab, Sb, Sbc, Sc and Scd, respectively.) The 
continuous line is the fitted result of the median value of each type (represented by the crosses). 


by using the method shown in this paper. However, these 
values are also in the range of the measurement error. 
Therefore, the method presented in Zhao, Peng & Wang 
(2004), which is more easily to be carried out, might be 
available and might be used to measure the thickness of 
the disks of nearly face-on galaxies for large samples. 
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Appendix A: Mathematical basis 

In order to study the effect of finite thickness of the disk on 
the dynamical properties of disk galaxies, it is convenient 
to adopt cylindrical coordinates (r, 9 , z) with the galactic 
center chosen at the origin of the coordinate system. The 
galactic plane is depicted by z = 0 and 9 represents the az¬ 
imuthal coordinate. For disk galaxies with zero-thickness, 


the self-gravitational potential <E> is governed by Poisson’s 
equation 

V 2 $ = 47rGcr(r, 9)5(z), (A.l) 

where G is the gravitational constant, V 2 denotes the 
Laplacian operator, S(z) is Dirac’s delta function that con¬ 
fines all the disk mass to the galactic plane and er(r, 9) rep¬ 
resents the surface density. In regions outside the galactic 
plane z ^ 0, equation (A.l) then becomes 

V 2 $(r, 9, z) = 0. (A.2) 

since there is no distribution of matter outside the disk 
plane. It is particularly convenient at this point to in¬ 
troduce the Laplace transform for the vertical coordinate 
z and the Fourier transform for the azimuthal angle 9 , 
namely 

/»oo 

A>(r,8,z)=e~ irne U 0 {r)e~ 0 ^dp. (A.3) 

J o 

Applying these transformations to equation (A.2) we then 
have 

x-^-{x dU0 ^ ) + (a; 2 - m 2 )Up{x) = 0, (A.4) 

ax ax 

where x = f3r, and the solution to equation (A.4) is the 
well known Bessel function of order nr 

Uf}{ x) = U((3r) = -J m (j3r). 


(A.5) 
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Using equation (A.5) in equation (A.3), we obtain References 


/*oo 

$(r, 0, z) = e~ im6 / [~J m ((3r)\ e~^dp. (A.6) 

Jo 

On the other hand, we may integrate Poisson’s equation 
(equation (A.l)) with respect to z to obtain the following 
expression for the surface density <r(r, 8) 


cr(r, 8) 


1 

UrG 


5$ 




o-}- 


(A.7) 


Substituting equation (A.6) in equation (A.l) and equa¬ 
tion (A.7) we can recast Poisson equation for disk galaxies 
with zero thickness in the form 


V 2 {/ 0 °° e~ ime {-J m {tJr)\ e~^dp) 
= f 0 °° e~ imff 2pj m (Pr)6(z)d{3, 


(A.8) 


V 2 [J m (pr)e- im0 -^~ z ' I] = -2pj m (Pr)e- tme S(z-z').(A.9) 


Appendix B: Solution of Poisson’s equation for an 
axisymmetric disk 

On the basis of the mathematical treatment just depicted 
together with the standard method of Green functions, it 
is straightforward to derive the appropriate Poisson equa¬ 
tion for the three-dimensional disk galaxies with finite 
thickness 


V 2 <p(r, z) = 4:TrGp(r, z ) 


(B.l) 


where the vertical distribution of matter is depicted by 
Parenako’s law, p(r, z) = ^a(r)e~ a ^ z ^ (see equation (3) in 
the text). To proceed, we apply the Bessel-Fourier trans¬ 
form to <j{t) to obtain 

/»oo 

a{r) = / pj 0 (pr)S(p)dp, (B.2) 

Jo 



rJo(pr)a(r)dr. 


(B.3) 


where S(/3) is the Bessel-Fourier transform for <r(r). 
Substituting equations (A.l) and (B.2) into equation 
(B.l), the three-dimensional Poisson equation can be 
rewritten as 
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V 2 </>(r, z) = 2nGa 


, , (B.4) 

■ fZo e ~ a \ Z '\dz' J 0 °° PJo(rl3)S(mz - z')df3. 


Compare equation (B.4) with equation (A.9) we then find 
the formal solution of the gravitational potential for a 
galactic disk with scale height H = 2/a 

r>OC 

ip(r, z) = — 27rGa / Ja(r/3)S(/3)F(a, (3, z)d/3, (B.5) 

Jo 


where 


F(a,P,z) 


1 

P 2 - a 2 


Pe- a 1*1 -ae-PW 


(B.6) 






